Spin textures in condensates with large dipole moments 



O 
(N 

1—5 



C/2 



I 

o 



> 

O 
(N 

l> 

O 
O 



X 



J. A. M. Huhtainaki,i'2 m. Takahashi/ T. P. Simula/^^ T. Mizushima/ and K. Machida^ 

^Department of Physics, Okayama University, Okayama 700-8530, Japan 

^Department of Applied Physics/COMP, Aalto University School of Science and Technology, 

P.O. Box 15100, FI-00076 AALTO, Finland and 

^School of Physics, Monash University, Victoria 3800, Australia 

We have solved numerically the ground states of a Bose-Einstein condensate in the presence of 
dipolar interparticle forces using a semiclassical approach. Our motivation is to model, in particular, 
the spontaneous spin textures emerging in quantum gases with large dipole moments, such as ^^Cr 
or Dy condensates, or ultracold gases consisting of polar molecules. For a pancake-shaped harmonic 
(optical) potential, we present the ground state phase diagram spanned by the strength of the 
nonlinear coupling and dipolar interactions. In an elongated harmonic potential, we observe a 
novel helical spin texture. The textures calculated according to the semiclassical model in the 
absence of external polarizing fields are predominantly analogous to previously reported results for 
a ferromagnetic F — 1 spinor Bose-Einstein condensate, suggesting that the spin textures arising 
from the dipolar forces are largely independent of the value of the quantum number F or the origin 
of the dipolar interactions. 



I. INTRODUCTION 

Long-range interparticle forces in a quantum system 
with a large coherence length is an intriguing combina- 
tion bound to exhibit a host of fascinating phenomena. 
Perhaps the most timely example of such a system is the 
gaseous atomic Bose-Einstein condensate (BEC) subject 
to magnetic dipole-dipole forces fll. 

The dipolar interaction potential, decreasing as r~^ 
in terms of the interparticle distance r, dominates on 
length scales determined by the coherence length. Other 
two-body interactions present in the system, such as in- 
duced dipolar forces (van der Waals), weaken typically 
much faster (r~^) and become negligible already over 
distances of an average interparticle separation. A fur- 
ther interesting aspect of the dipole-dipole interaction is 
its anisotropy enriching the already diverse finite-size ef- 
fects in trapped ultracold atomic gases. The magnetic 
dipolar interaction in condensates has been predicted to 
give rise to phenomena ranging from spin textures and 
spontaneous mass currents [2tJ| to roton minimum in the 
excitation spectrum [5|, |6[ , linking the field into the study 
of liquid He II. 

The realization of ^^Cr condensates has provided 
means of probing dipolar effects experimentally due 
to the exceptionally large magnetic moments of the 
atoms [3| ■ The ground states of a chromium condensate 
have been studied extensively |8l-[l3|. Anisotropic de- 
formation of an expanding chromium condensate due to 
dipolar forces has been observed llM j ^^id dipole- induced 
spin relaxation in an initially polarized ^^Cr has been 
linked to the famous Einstein-de Haas effect in ferro- 
magnets 0, llJI- Also, collapse and subsequent d-wave 
symmetric explosion of dipolar condensates have been 
recently studied in the case of ^^Cr both experimentally 
and theoretically [lal • Chromium condensates have been 
recently produced through optical methods |la |. 

The strength of the magnetic dipolar interaction is de- 
termined by the atomic magnetic moment ^im through 



the coupling constant 5^ = ^q^^j/Att, where /io is the 
permeability of vacuum. For example, for alkali con- 
densates with total angular momentum quantum number 
_F = 1, the magnetic moment is given by ^m = fJ'Bgp, 
where /i^ is the Bohr magneton and gp ~ 1/2 the 
Lande (7-factor. Such systems are subject to weak dipo- 
lar interactions, e.g., g'^j/g' ^ 10^'^ for ^^Rb, where 
g' = Airh^ [ag + 2a2) /3m is the mean-field density- 
density coupling constant. Here ag and 02 are the s- 
wave scattering lengths in the channels with total spin 
and 2, and m is the atomic mass. Nevertheless, dipolar 
effects have been predicted to be observable in i^ = 1 
alkali BECs even in the presence of a magnetic field [l3l , 
which was recently confirmed experimentally based on 
time-evolution study of a helical spin texture [I3] ■ It has 
also been proposed that spin echo in spinor BECs could 
be utilized in revealing dipole-dipole interactions jl9| . 

The spontaneous occurrence of novel ground-state spin 
textures in the absence of external magnetic fields re- 
quires typically stronger dipolar interactions, g'j^/g' ~ 
10~^-10~^. Hence, the ^^Cr condensates consisting of 
particles with magnetic moments of 6 ^b , as opposed to 
maximal magnetic moments of 1 /is in alkali gases, seem 
more favorable for observing such effects. Moreover, the 
rare-earth-metal element Er, with a magnetic moment 
oi7 hb, has been cooled down to ^K temperatures [20| . 
Also, recent developments in trapping and cooling of Dy 
with the largest atomic magnetic moment of 10 iib yields 
a promising candidate for observing the predicted spin 
textures 2l|. Developments in the study of ultracold 



polar molecules provides means of investigating dipolar 
effects with large electric moments |22h24| . 

The study of alkali condensates based on a quantum 
mechanical mean-field treatment predict spin textures 
with the smallest possible value for the total angular mo- 
mentum quantum number with internal degrees of free- 
dom, namely _F = 1. It is worthwhile to approach an 
analogous problem from the other extreme limit by treat- 
ing the magnetic moments of the gas classically Q, [23 . 



By comparing the results predicted by the two models, 
one may expect that if the predictions agree, they could 
be of universal character for all dipolar condensates and 
independent of the particular value of the quantum num- 
ber F. In general, the quantum mechanical order param- 
eter has 2F + 1 components and the short-range interac- 
tion term contains F+l independent coupling constants. 
Hence, it would be very cumbersome to treat each value 
of F separately with the complexity of the problem in- 
creasing along with F. 

Our semiclassical model is briefly described in Sec. II. 
The order parameter is written in an alternative form 
compared to previous studies @, UM i'^ order to simplify 
analysis and to increase numerical efficiency. The main 
results are explained in Sec. Ill: The ground states of 
the system in harmonic traps of various geometries are 
described and the collapse of the spin vortex state is an- 
alyzed briefly. The novel spin helix state is introduced 
before concluding remarks of Sec. IV. 



II. MODEL 

In this Section, we construct a phenomenological 
mean-field model describing a trapped Bose-Einstein con- 
densate with local as well as nonlocal interparticle inter- 
actions. The model is equivalent to the semiclassical ap- 
proach previously studied in [2, [2^ , with the exception 
that now the order parameter field is written in a carte- 
sian basis yielding a set of three Gross-Pitaevskii type of 
equations leading to more efficient numerics. 

The local interaction is assumed to be of the stan- 
dard s-wave form, with coupling constant g' = Airh'^a/'m, 
where a is the s-wave scattering length. Henceforth, we 
will refer to its dimensionless form g = AirNa/ar, ex- 
pressed in natural trapping units: hiOr is the unit of en- 
ergy with uJr being the radial trapping frequency of the 
confining harmonic potential, and the radial harmonic 
oscillator length a^ = \fhJmiOT. is the unit of distance. 
The number of confined atoms is denoted by N . The 
nonlocal interaction is the anisotropic dipole-dipole inter- 
action with the dimensionless coupling constant 5^ which 
is quantified in relation to g throughout the article. Re- 
gardless of whether the origin of the dipolar interactions 
is considered to be magnetic or electric, we adopt no- 
tation and terminology as if it were of the former. For 
the trapped system to be stable, we find that for strong 
enough contact interaction, the value of g^ should not ex- 
ceed '-^ gl'^t a value very close to the number calculated 
in the F = \ case [3]. 

The order parameter is taken to be a three-component 
real- valued vector '\\) = (V'a:,'(/'yj V'z)- It is straightfor- 
ward to show that all line defects in such order-parameter 
space are topologically unstable, because any closed 
curve on a sphere can be continuously transformed into 
a point which corresponds to a spin polarized state [26j. 
Nevertheless, the energetically stable states can have 
nontrivial spin textures. 



In the present model, the particle density is assumed 
to be related to the order parameter through n(r) = 
^j,V'fc(r) and is normalized to unity, / drn(r) — 1. Let 
us make the assumption that the system is ferromagnetic, 
and hence we may require that all spins are pointing into 
the same direction within a small enough region of space. 
With this simplification, the magnitude of the magneti- 
zation density is related to the particle density through 
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(1) 



where fiM is the magnetic moment of a single particle, 
and Mk are the components of magnetization. By squar- 
ing, we obtain T^k^ki^) = Mm Efc V'fe(r)»^(i"), which is 
satisfied if we define 



Mfc(r) = ^iMil^k{r)\/n{r), 



(2) 



relating the magnetization density to the order parame- 
ter. In the following, we omit writing the constant fj,M 
explicitly and assume it to be included in the coupling 
constant gd- 

The energy functional i?tot['0a;j'0i/j'0z] can thus be 
written as 

Etot = JY. "Pkhi^kdv +11 n^dv + Y f [ ^('■' r')rfrdr', 

k 

(3) 
where h = — -^V^ + Vtraplr) is the single-particle Hamilto- 
nian and Vtrap = | {x^ + y'^ + A^z^) is the external trap- 
ping potential expressed in natural trapping units. For 
now we will omit external rotation and mass currents in 
the system, and hence the kinetic energy is merely due to 
quantum pressure. The second term in Eq. ^ describes 
the local mean-field s-wave interaction with the coupling 
constant (7, and the final term the nonlocal dipole-dipole 
interaction with 

D(r, r') = [M(r) • M(r') - 3 (M(r) • e^) (M(r') ■ e^)] /i?^ 

(4) 
where R = r — r' is the relative coordinate and e^ the 
unit vector along it. 

Stationary states of the condensate are obtained by 
differentiating the energy functional with respect to the 
components of the order parameter with the particle 
number constraint taken into account through a La- 
grange multiplier /i. Differentiation with respect to ^pj 
results in a set of three Gross-Pitaevskii equations 



hipj + gntpj -\ 



E, Muh 



ipj + y/nlj 



= l^i'j 



(5) 



Here the functions Ij are defined by 



/.(r)^/ 



Af,(r')-3e],^A/,(r')e 



/R^'dv', (6) 



with e'^ being the /th component of e^. These integrals 
may be further broken into convolutions. By applying 



the convolution theorem, we obtain 



Ij{v)^F- 



^J-[MKr)]/y(k) 



(7) 



where T stands for Fourier transform and 
//j(k) = —^ {5ij — 3kikj/k^) is the Fourier trans- 
form of //j(r) — [5ij — iriTj/r"^) jr" . The Fourier 
transforms are efficiently evaluated by using Fast Fourier 
Transform. 

From the general form of the GP equations, Eq. ([S]), 
it is possible to conclude that the spin-polarized texture 
is not a stationary state in a confined three-dimensional 
system in the absence of external polarizing fields when 
ga ^ 0. Namely, Eq. dS]) is of the form Ajk^pk = bj, where 
bj — ^-y/nlj ^ 0, in general. In the spin-polarized state, 
we may choose, say, the z-axis along the polarization, 
whence ip^ — ipy = yielding bx — by = from the 
general form above. When gd ^ 0, this can be satisfied 



in regions of non- vanishing density only if I^ = ly 



0. 



Hence, the bracketed expression in Eq. ([71) must vanish 
identically. As jzx and j^y are non- vanishing in any finite 
volume Sk, continuity of f/'^ implies M^{x) = 0, which 
is a contradiction. Such conclusion can also be drawn 
from the quantum mechanical model by following similar 
arguments. 

Apart from the quantum pressure term in Eq. ([3]), 
the spin model described above can be viewed as re- 
sulting from a classical energy functional. However, 
the present model can also be argued from the quan- 
tum mechanical spin-_F model constrained within the 
ferromagnetic manifold [27| . With maximally aligned 
spins, the order parameter at a fixed point r is of 
the form -0 = ^e^e^-^F.<^/fi(,-^FyP/'i(,-'F.i/'^\z) = 
^^■i{e-F',)^-iF.a/h^-iFyfi/h^^^^ where F^ are the hyper- 
fine spin operators, and Fz\z) — hF\z). The order param- 
eter of the classical spin model is obtained if we neglect 
the phase factor e*'^*"^^-' and replace the quantum me- 
chanical rotation operators by the classical equivalents 
and the eigenstate \z) by the unit vector pointing along 
the z-axis. Such substitution should be valid when quan- 
tum fiuctuations of the spin operator F become negligi- 
ble. The relative fiuctuations in the state \z) are given 
by ((F - hFf)/{Y'^) = 1/{F + 1), which vanish in the 
limit of large F . Possible mass currents arising from lo- 
cal spin-gauge symmetry are neglected when e^^^~^^' is 
set to unity, and the kinetic energy reduces merely to the 
quantum pressure term in Eq. ([3]). 



III. RESULTS 

We have solved the ground states of the system with 
various values of the coupling constants g and gd, and the 
aspect ratio A. Special emphasis is given to the pancake- 
and cigar-shaped systems, for which we choose A = 10 
and A = 0.10-0.50, respectively. 




FIG. 1: (Color online) Magnetization M(r) (arrows) and den- 
sity n(r) (color) of the (a) flare and (b) spin vortex states for 
g = 100, gd/g = 0.15 in a trap with aspect ratio A = 10. Both 
quantities are shown in the z = plane, the density being 
nearly gaussian in the axial direction and M^ small. For the 
chosen parameter values, the flare state in (a) is the energeti- 
cally favored conflguration. Panel (c) illustrates a spin vortex 
state with opposite spin winding compared to (b). Such state 
is not energetically favorable for the parameter values consid- 
ered in this work. Each panel has dimensions Sur x Sur- 



A. Ground states in the pancake-shaped limit 

Let us first consider the case of a cylindrically sym- 
metric harmonic trap with strong, A = 10, confinement 
in the axial (z) direction. In the presence of dipolar in- 
teractions, gd > 0, the magnetic moments tend to lie 
predominantly in the plane perpendicular to the axial 
direction in order for the system to minimize dipolar in- 
teraction energy. 

For small enough value of gd/g, the spin texture has 
typically the flare structure which has been studied pre- 
viously using the semiclassical approach as well as the 
quantum mechanical mean-field model in the F == 1 
case @tJ|- Such state is illustrated in Fig. [IJa) for, 
g = 100 and gd/g — 0.15. The arrows denote the lo- 
cal direction of magnetization M(r), whereas the color 
refers to the particle density n{r). The repulsive interac- 
tion between parallel spins separated by a vector perpen- 
dicular to the spin vectors causes the magnetization to 
deviate from the polarized texture. The structure may 
also be thought of as resulting from the presence of two 
spin vortices located at the periphery of the cloud. The 
spin texture is flare-like also in the x-z -plane (y = 0) 
due to finite M^ , which is in accordance with the picture 
that a single toroidal spin vortex encircles the cloud. In 
the flare state, the magnetization has even parity. 

When the strength of dipolar interactions is increased, 
the ground state undergoes a second order phase tran- 
sition into a state hosting a single spin vortex which is 
iUustrated in Fig.IHb) for g = 100 and gd/g = 0.15. The 
density is typically suppressed at the core of the vortex. 
The spin vortex state has also been studied previously 
within the semiclassical as well as the F = 1 case 12tJ1- 
Analogously to the flare state, the presence of the spin 
vortex results in a texture which favors dipolar interac- 
tions by reducing the repulsive interactions of parallel 
spins separated by a vector perpendicular to their mag- 
netization. For example, close to the phase transition 
line in Fig. H with g = 1000 and gd/g = 0.05, the differ- 



ences in the kinetic, potential, contact interaction, and 
dipolar energies of the flare and the spin vortex states 
are AE^^in = -0.f6, A^pot = 0.077, A^ni = -0.11, and 
Ai^dip — 0.28, respectively, leading to a gain in total 
energy of Aii'tot = 0.077 in units of hujr per particle. 

Figure [T]Jc) illustrates a spin vortex state with oppo- 
site spin winding. For this texture, the angle between lo- 
cal magnetization and the cc-axis decreases as the vortex 
core is circled around in the counterclockwise direction, 
whereas for the state in Fig. [Tfb), the angle increases. 
Such state is found only as an excited solution in the 
present work. In a larger dipolar system, one can con- 
struct energetically low-lying spin vortex lattices by ar- 
ranging the vortices presented in Figs.HJb) andfljc), and 
their negative counterparts (M — >■ — M) in an alternat- 
ing square lattice. Both spin vortices presented in Fig. [1] 
have odd parity. 
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FIG. 2: Ground-state phase diagram of a dipolar condensate 
in a harmonic trap with aspect ratio A = 10. The effective 
contact interaction coupling constant g is represented in log- 
arithmic scale on the vertical axis. The horizontal axis, mea- 
suring the strength of dipolar interactions through the ratio 
Qd/g, has linear scale. The phase diagram is divided into 
three regions: flare (Fig. Ula)), spin vortex (Fig. [Ub)), and 
the region where the spin vortex becomes unstable against 
collapse. 

The ground-state phase diagram in the {g,gd/g)- 
parameter plane is shown in Fig. [2] The axes of the plane 
are chosen such that the abscissa is proportional to the 
particle number N and the ordinate is independent of 
N and proportional to the (bare) dipolar coupling con- 
stant g^. The spin vortex state is energetically favored 
for strong contact and dipolar interactions. The flare 
state dominates the phase diagram in the limit of weak 
contact interaction, 5 < 50, regardless of the strength 
of dipolar interactions. The phase transition point from 
flare to spin vortex state depends strongly on the value 
of g. The spin vortex state becomes unstable towards 
collapse beyond the critical value of gd/g ~ 0.25 which 



depends only weakly on the value of g for g > 500. 

With strong enough contact interaction {g > 1000) and 
weak dipolar interaction, the cores of the spin vortices are 
filled with particles whose magnetic moments are point- 
ing in the axial direction. The condensate gains trapping 
as well as contact interaction energy by filling the vortex 
core. The dashed line in the phase diagram of Fig. [5] sep- 
arates the states with filled cores from states with empty 
cores. On the line / n{r — 0, z)dz/maxj n{r, z)dz — 0.01 
with r = y'x^ + y^, whereas the ratio is close to unity in 
the upper left corner. The finite axial magnetic moment 
due to the filled core breaks the inversion symmetry of 
the state. Instead, the components of M have the follow- 
ing symmetry: PzM^^y = -M^^y, P^M^ = M^, where P^ 
inverts the sign of the x- and y-coordinates keeping z 
intact. 

For large enough g, the flare state develops continu- 
ously into a state with two spin vortices which have fer- 
romagnetic cores as the strength of dipolar interactions is 
increased. The magnetic moments of the cores are point- 
ing either into the same or opposite directions, the two 
states being nearly degenerate irrespective of the rela- 
tive orientation. States hosting multiple spin vortices are 
found to be energetically unfavorable compared to single 
spin vortex states for the parameter values considered in 
this work. 

The radial size of the spin vortex state diminishes sig- 
niflcantly as the strength of dipolar interactions is in- 
creased. This suggests that the reason why the system 
becomes unstable at some critical value of gd/g could 
be due to inward collapse of the condensate. Local and 
global collapse of a dipolar condensate has been recently 
studied numerically [2S] . 

In order to understand why the spin vortex solution 
ceases to exist above the critical point, it is instructive 
to study scaling transformations of the form 

f^{T)Mr, z) = c"(T)V'fc([l + T]r, [1 + aT]z), (8) 

where (r, z) are the cylindrical coordinates, t is the scal- 
ing parameter, a determines the ratio between axial and 
radial scaling, and c"' (r) is chosen to ensure particle num- 
ber conservation. 

Close to the critical point of collapse, the spin vor- 
tex state in a pancake-shaped trap is the ground state of 
the system, and hence lies in a minimum of the energy 
functional. Under transformations of the form given in 
Eq. ([U , the total energy becomes a function of the scal- 
ing parameter r, E^^^{t) = Etot[T''{T)ipk{r, z)]. Devia- 
tion from the ground state always leads to increase in 
energy, and hence the second derivative of the total en- 
ergy with respect to any one-parameter transformation 
must be positive, 9^i?^t(r)|^ > 0. The existence of a 
transformation for which this quantity vanishes indicates 
that the state becomes unstable against such variation. 

Figure [3] shows the value of m.in„{^^E^^^.{T)\^ } as 
a function of gd/g S [0.10,0.30] scaled by the value at 
gd/g = 0.10. The solid curve is for g = W^, dashed for 



g — 10^, and dash-dotted for g = 10^. The curves are 
extrapolated (dotted lines) using the last few points to 
obtain an estimate for the critical value for which the 
minimum in the energy functional vanishes. The criti- 
cal values are gd/g = (0.240,0.243,0.274), respectively. 
The inset in Fig. [3] depicts the value of cr for which 
9^£'^j(t) is minimized for each gd/g, the horizon- 

tal axis being the same as in the main graph. In the 
vicinity of the critical point for g < 10'^, ct > 0, showing 
that the collapsing cloud shrinks both in radial and axial 
directions. 
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FIG. 3: (Color online) Second derivative of the total energy 
with respect to a scaling transformation of the form given 
in Eq. ((8]) as a function of gdl 9 shown in units of the cor- 
responding quantity at gdl g = O.fO. The curves correspond 
to the parameter values g — 10* (solid), g = 10'^ (dashed), 
and g = 10'^ (dash-dotted). The value of gdj g for which the 
second derivative vanishes indicates the critical strength be- 
yond which the spin vortex state becomes unstable against 
collapse. The inset shows the ratio of axial and radial scaling 
for which the minimal value of the bracketed expression in 
the main figure is obtained. 



Based on a spinor F — \ study, a critical value of 
9dl 9 ~ 0.24 has been previously reported for the exis- 
tence of the spin vortex state [3'] , where the parameters 
are chosen such that g w 7000 in the present study. As 
the current model is expected to be accurate, apart from 
possible mass currents, in the limit of large (classical) 
magnetic moments, this agreement suggests that the crit- 
ical value for the collapse is universal and independent of 
F. 



with radial frequency w^ = 27r x 100 Hz. Regardless of 
the aspect ratio A, the solutions are found to exist only 
within the interval gdl 9 G [0,0.235], agreeing|with the 
result previously reported in the F — \ study '3'|. 

The spins tend to lie predominantly along the axial 
direction for finite but sufficiently weak dipolar interac- 
tions. This ground state resembles the flare state in the 
pancake-shaped limit, and it has been discussed previ- 
ously both in i^ = 1 condensates as well as using the 
semiclassical model [2,141. Figures |4l^a)-(c) illustrate the 
spin textures in three radial cross-sections of the conden- 
sate. Here A = 0.20, gdj g = 0.030, and the cross-sections 
are taken at z = — 12ai.,0, and 12 a^, respectively. The 
color depicts the z-component of magnetization, M^iv)^ 
and the color bar is scaled with respect to the maximum 
magnetization, max|M(r)|, in the corresponding state. 
The arrows show the texture projected onto the x-y - 
plane (henceforth referred to as the planar texture), with 
the length of the arrows scaled within each panel sepa- 
rately, except in Fig. Hfb) , for which M^ — My = by 
symmetry. The stability range of the flare state depends 
strongly on the aspect ratio: For example, with A = 0.50, 
the flare state is the ground state for < gdl 9 ^ 0.01, 
with A = 0.20 for < gdl g < 0.08, whereas with A = 0.10 
the flare state dominates the entire stability window. 

Typically, for stronger dipolar interactions, a spin vor- 
tex texture appears in the central region of the con- 
densate, illustrated in Figs. |ll[d)-(f) for A = 0.20 and 
gdig = 0.080. As shown in Fig. life), the planar tex- 
ture near the center resembles the spin vortex texture 
in the pancake-shaped geometry discussed before, see 
Fig. [IJb). As one moves further away from the center, 
the planar texture deforms continuously towards the flare 
texture discussed in the previous paragraph, as depicted 
in Figs, md) and|3Jf). Although the magnetization in 
both the flare and the spin vortex states has the same 



symmetry, PzMx,. 



-M^^y, P,M, 



Mz, the phase 



transition is sharp, as illustrated below in Fig. [5jd). 

In order to characterize the spin vortex state more pre- 
cisely, we deflne the following quantities: The axial col- 
umn density reads 



nz{z) = / n{r)dxdy. 



(9) 



This measures the number of atoms per unit length in the 
axial direction and is normalized to unity. The average 
twisting angle is given by 



a{z) — / arccos 



>-xy 



■M,,(r) 



M.y{r) 



(10) 



B. Ground states in the cigar-shaped limit 

Let us now consider solutions to Eq. ([5]) in an elon- 
gated trapping geometry with A e [0.10,0.50]. For def- 
initeness, we will Hx g — 10* which corresponds to a 
number of N ^ 1.5 x 10^ *^Rb atoms in a harmonic trap 



where r^y — {xe^ -V yey) I \/ x"^ +2/^ and the averaging 
is taken over vectors Ma;^ whose length exceeds 1 % 
of the maximum of the planar magnetization M,j.y = 

A/2 -f M^ . This quantity characterizes the twisting 

of the magnetization in plane, yielding zero (or tt) for 
the fiare-like textures. Figs. IHa)-(c), and 7r/2 for the 
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FIG. 4: (Color online) Spin textures in the flare, (a)-(c), 
spin vortex, (d)-(f), and spin helix, (g)-(i), states in a cigar- 
shaped trapping geometry with aspect ratio A = 0.20 and 
dipolar interaction strengths Qd/g = 0.030, 0.080, and 0.030, 
respectively. The arrows illustrate the magnetization within 
a given radial cross-section projected onto the x-y -plane, 
whereas the color refers to the axial magnetization M^ nor- 
malized with respect to maximal magnetization within each 
state separately. Each panel has dimensions 12 a^ x 12 ar. 



spin vortex texture. Fig. EJe). The average tviristing an- 
gle is essentially independent of the radial distribution of 
the density. Finally, we define the average tilting angle 
through 



(3{z) — / arctan 



M.y{r) 



M,{r) 



(11) 



vi^here the averaging is evaluated as above. The tilting 
angle is related to the pitch of the helical streamlines 
obtained by following the local direction of magnetization 
in the spin vortex state, c.f. Ref. [3]. 

Figure[5l^a) shows the twisting a{z) in the flare (solid) 
and spin vortex state (dashed and dash-dotted) for the 
dipolar interaction strengths gd/g = (0.075,0.15,0.235), 
respectively. The flare state is chosen from the neigh- 
borhood of the transition point to a spin vortex state. 
However, a(z) remains nearly zero (or tt) over the whole 
length of the cloud. Small deviation from zero shows 
that the flare state has even parity only approximately. 
As the dipolar interaction strength is increased, a spin 
vortex enters the system. Hence, the twisting angle de- 
creases continuously from tt to along the length of the 
condensate. A plateau of a{z) « tt/2 forms in the cen- 
tral region of the system for strong dipolar interactions, 
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FIG. 5: (Color online) The twisting angle a(z) (a), tilt- 
ing angle I3{z) (b), and the axial column density nz{z) (c) 
for three values of the dipolar interaction strength gd/g = 
(0.075, 0.15, 0.235) shown with solid, dashed, and dash-dotted 
lines, respectively. The total magnetization Mtot, which is di- 
rected along the z-axis by symmetry, is shown in (d) for the 
aspect ratios A = 0.10 (dashed), 0.20 (solid), and 0.50 (dash- 
dotted). The dots in (d) refer to the values of gd/g used in 
(a), (b), and (c). 



gd/g ^ 0.15. The width of the plateau decreases for in- 
creasing gd/g due to shrinking of the cloud. 

The tilting angle /3{z) is shown in Fig.ISJb) for the same 
parameter values as in Fig. [Sfa). It remains relatively 
small in the flare state and experiences a sudden increase 
at the center of the system when the ground state hosts a 
spin vortex. For very strong dipolar interactions, I3(z) w 
7r/2 in the central region, slightly even exceeding 7r/2 
due to interaction with the axial magnetization of the 
core region. The small lobes in /3{z) close to the top and 
bottom of the cloud are remainders of the flare state. 

The strength of dipolar interactions affects the spatial 
density profile of the spin vortex state significantly. The 
axial column density nz{z) is shown in Fig. [SJc) for the 
parameter values used in 5(a) and 5(b). Not only does 
the system shrink in the radial, but also the axial di- 
rection with increasing gd/g- Also, the column density 
appears to be slightly bimodal for strong enough dipo- 
lar interactions: the column density is enhanced in the 
central region of the condensate where the spin vortex 
lies in order for the system to gain dipolar energy. The 
width of the plateau in Fig. [5ja) due to the presence of 
the spin vortex matches the size of the central profile 
in the bimodal density distribution. The bimodality ap- 
pears more vividly in elongated systems with A < 0.20. 
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FIG. 6: (Color online) (a) Mx (solid) and Aly (dashed) on the ^-axis as a function of z in the spin helix state for A = 0.10 
and Qd/g = 0.20. The dash-dotted curve in the inset shows M^ along the y-axis in the range [—5,5]ar. The vertical axis 
spans the interval [—0.0028, 0.0028] N^m /ai- both in the main figure and the inset, (b) The column density n{x, z) = j n{r)dy 
illustrating density oscillations characteristic of the spin helix state for strong dipolar interactions. The parameters are as in 
(a), and the field of view is 8 a^ x 60 a^. (c) Wave vector of the spin helix state for Qd/g ~ (0.080, 0.20, 0.23) shown with solid, 
dashed and dash-dotted curves, respectively. The peaks at the ends of the cloud are finite-size effects. 



In the extreme limit of gd/g — 0.235, the density n(r) in 
the spin vortex state is significantly reduced close to the 
center of the trap where r « 0. 

Figure [5jd) depicts the total magnetization Mtot = 
J M(r)(ir in the flare and spin vortex states as a func- 
tion of gd/g- Due to symmetry, the total magnetiza- 
tion is along the axial direction. The dashed, solid 
and dash-dotted lines correspond to the aspect ratios 
A = (0.10,0.20,0.50), respectively. The sudden drops 
in total magnetization indicate the phase boundary be- 
tween the flare and spin vortex states. Whereas the flare 
state dominates the entire stability window for very elon- 
gated trapping geometry (A — 0.10), the ground state in 
a prolate system (A = 0.50) hosts a spin vortex already 
with gd/g = 0.01. The phase transition points for dif- 
ferent values of A agree qualitatively with the analogous 
results in the F — I study [3] . 

It is reasonable to expect that there exist also station- 
ary states with the opposite symmetry compared to the 
flare and spin vortex states, i.e., PzMx,y{r) — M^^y^r), 
PzMz{v) — -~Mz{v). There indeed exist low-energy so- 
lutions to Eq. ([5]) with such symmetry, to which we refer 
to as spin helices. The spin helix state is found, e.g. 
with A = (0.10,0.20,0.50) in the entire stability interval 
< gd/g ^ 0.235 of the system. This state resembles 
closely the state studied in Ref. [l^, where the helical 
spin texture is created by using a transient magnetic field 
gradient. According to our simulations, the stationary 
spin helix state exists also in a ferromagnetic F = 1 sys- 
tem with dipolar interactions, which will be studied in 
more detail elsewhere. Dynamical instability of a similar 
structure in the absence of dipolar interactions has been 
studied recently |29j . 

The spin helix texture is illustrated for A — 0.20 and 
9d/g = 0.030 in Fig. |llJg)-[Hi), where the radial cross- 
sections are taken dX z = (—2.5 a^, 0, 2.5 a^), respectively. 
On the 2-axis, the magnetization lies in the x-y -plane, 



as a consequence of the antisymmetry of M^. Further 
away from the z-axis and perpendicular to the magneti- 
zation on the axis, M^ becomes the dominant component. 
The whole planar texture rotates about the z-axis as a 
function of the z-coordinate, traversing typically through 
several cycles along the length of the condensate. 

Energetically, the spin helix state appears to be fa- 
vored by strong dipolar interactions and not too elon- 
gated geometries. For example, with A — 0.50, the he- 
lix becomes energetically favorable compared to the spin 
vortex state between 0.050 < gd/g < 0.10. It is challeng- 
ing to pinpoint the exact location of the phase transition 
point due to near degeneracy of the two states. Near 
the critical value of gd/g — 0.235, the difference in the 
total energy between the spin vortex and helix states is 
A£'tot ~ 0.1 hujr per particle in favor of the helix, which 
is roughly 2% of the total energy. For A = 0.20, the helix 
state appears to be the minimal energy texture only for 
gd/g ^ 0.20, whereas for A = 0.10, the flare state lies 
2%-7% lower in energy over the entire stability range of 
the system. 

The number of cycles in the helix texture increases 
as the condensate is elongated, and thus for the sake of 
clarity we illustrate it as an excited state for A = 0.10 and 
gd/g = 0.20 in Figs.inija)-(c). In Fig. IHlJa), the solid and 
dashed curves show the components M^ and My along 
the z-axis, respectively, and the Mz component along the 
y-axis is shown in the inset by the dash-dotted curve. 

The spin helix can be thought of as two elongated 
stripes, polarized along the z-axis in the opposite di- 
rections, intertwined around one another. The helical 
texture on the z-axis arises due to continuous twisting 
of the magnetization via the x-y -plane. In the F — 1 
case, quantized spin vortices of opposite winding pene- 
trate through the axially polarized ferromagnetic stripes, 
forming an intertwined spin vortex pair. Intertwining of 
two mass vortices has been previously studied in relation 
to the splitting of a doubly quantized vortex in a scalar 
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condensate [30, ^V\ . 

Figure ^h) illustrates the column density n{x, z) — 
J n{r)dy for the same state as in Fig. |6l[a), red denot- 
ing area of high and blue of vanishing particle density. 
With strong dipolar interactions density oscillation ap- 
pear spontaneously due to the helix spin texture: for a 
fixed point in the x~y -plane close to the surface region 
of the cloud, the axial magnetization Mz is an oscillating 
function of z. The particle density is suppressed in the 
vicinity of the nodes of M^ and enhanced at the anti- 
nodes due to dipolar interactions. 

As a measure of the pitch of the spin helix, we define 
the angle 



(z) — arctan 



My{0,0,z) 



M,(0,0,z) 



(12) 



The derivative dO/dz yields the wave vector of the he- 
lix, which is plotted in Fig. ^c) for A = 0.10 and 
gjg = (0.080,0.20,0.23) with the solid, dashed, and 
dash-dotted curves, respectively. The wave vector tends 
to increase for stronger dipolar interactions, which is rea- 
sonable because the dipolar coherence length decreases as 

— 1/2 

Cd oc (7^ . The peaks in dO /dz at the top and bottom 
of the cloud are finite-size effects: The texture may ad- 
just freely into an energetically favorable configuration at 
the edge as one of the boundary conditions due to con- 
tinuity of the order parameter is liberated. Oscillations 
penetrate along the whole length of the condensate for 
9d/ g = 0.23. These oscillations enhance rapidly as the 
strength of dipolar interactions is increased even further. 
The number of cycles in the spin helix state decreases 
as the aspect ratio A — tOz/uir is increased, until in spher- 
ical geometry, the direction of the spin on the z-axis 
twists only through half a cycle along the length of the 
system, c.f. Fig.lHl^a). Interestingly, the spin vortex state, 
for which the magnetization has the opposite symmetry 
with respect to inversion about the z-axis, reduces to the 
spin helix state, rotated by (±)7r/2 about Gz x Mh, where 
Mh is the magnetization at the trap center in the helix 
state. 



IV. SUMMARY AND CONCLUSIONS 

We have studied spin textures arising from dipolar in- 
teractions in gaseous Bose-Einstein condensates of parti- 
cles with large permanent dipole moments. The theory 
is based on a semiclassical model treating the dipole mo- 
ments of the bosons classically. 

The observed spin textures in clouds confined in har- 
monic trapping potentials agree qualitatively with previ- 
ously reported results for an F = 1 system [3, Q] , such as 
^^Rb. Moreover, the ground-state phase transition points 



with respect to the strength of dipolar interactions seem 
to agree roughly both with weak and tight axial trapping 
frequency. The qualitative agreement in the observations 
drawn from the two models suggests that similar textures 
and phase diagrams are to be expected also for ferromag- 
netic systems with _F > 1 and for condensates consisting 
of electric dipoles. 

A major difference between our semiclassical model 
and the quantum mechanical model describing a mag- 
netic system lies in the existence of mass currents: In the 
quantum mechanical case, phase gradients of the com- 
ponents of the order parameter emerge spontaneously 
possibly giving rise to mass currents which are absent 
in the semiclassical treatment, c.f. the last paragraph of 
Sec. II. The main effect from taking spontaneous mass 
currents into account regarding the present results would 
be that the phase transition line between the flare and 
the spin vortex states in Fig. [2] is shifted to the right due 
to increased kinetic energy of the spin vortex state. Spin 
dynamics of ferromagnetic condensates has been studied 
recently in thelong-wavelength limit using a hydrody- 
namic model [32!|. 

In addition to solving the ground states of the sys- 
tem for various parameter values, we investigate the spin 
vortex state in the extreme limit of the dipolar inter- 
action strength gal g ~ 1/4. For larger dipolar interac- 
tions, the state becomes unstable against collapse of the 
cloud due to strong attractive forces overwhelming the 
quantum pressure term and repulsive interparticle inter- 
actions. The estimated point of instability agrees well 
with the value observed in the F = 1 study Q. 

In the limit of tight axial confinement, two ground 
states are observed, namely, the flare and the spin vortex 
states. For prolate geometries, an additional spin helix 
texture appears as a low-energy stationary state. The 
helix is the ground state of the system only in slightly 
prolate condensates and for strong dipolar interactions. 
This state is most likely related to the S state reported 
in [3], and is especially interesting in relation to the ex- 
perimental observation of dipolar effects in ^^Rb utilizing 
a similar spin texture [l8|. The magnetization pattern 
of the spin helix gives rise to spontaneous density oscil- 
lations in the stationary state for strong dipolar interac- 
tions. As in the case of the spin vortex and flare textures, 
the helix state ceases to exist for gd/g > 1/4. 
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